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I. Research Carried Out. 

The following report, which is being prepared for submission to the 
Astrophysics Journal, is a full description of research completed under 
tills grant. 



Observations of Comet Kohoutek (1973f) with 
a Ground-Based Fabry-Perot Spectrometer 

D, Huppler, R. J. Reynolds, F. L. Roesler, F. Scherb and 0, Trauger 
Physics Department, University of Wisconsin, Madison 

ABSTRACT 

Observations of Ha, and [01] emission Tines from comet Kohoutek 

were made during a two month period around perihelion. Analyses of Ha line 
profiles and line intensities indicate that the mean outflow velocity of 
the hydrogen atoms was 7.8 ± 0,2 km s ^ and that the hydrogen atom production 
rate varied from about 1.0 x 10^® s'^ to about 3.5 x 10^® for comet-sun 

distances between 1 AU and 0.4 AU, respectively. The identification of an 
HgO'^ emission feature in certain Ha scans indicates that the ions were 
moving in a tailward direction with a velocity of 20 to 40 km s"'' with respect 
to the_com‘et nucleus. An upper limi t of 1 part in 100 was found for the D/H 
ratio in the cometary atomic hydrogen cloud. 

Subject heading: comets 

I. INTRODUCTION 

Between 1973 December 1 and 1974 February, 2, optical emission lines from 
the gas cloud surrounding Comet Kohoutek were observed using a double Fabry- 
Perot etalon spectrometer at Kitt Peak National Observatory. The spectrometer 
had a resolving power of 40,000, corresponding to a velocity resolution of 
about 7.5 km sec T With this resolution it was possible to use the comet- 
earth relative velocity to resolve faint cometary Ha X6563, [01] X6300 and 
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other emission lines from geo.coronal and airglow emissions and to study the , 
.cometary line profiles in order to obtain information about the composition, 
effective temperatures, outflow velocities, and production rates of atoms 
and ions in the cometary envelope. . 

II. INSTRUMENTATION 

The spectrometer was coupled to the McMath Solar Telescope by focussing 
the primary image of the sky, which has a scale of about 1 arc min per 25 mm, 
directly onto the 150 mm diameter etalon of the Fabry-Perot spectrometer. 

Masks could readily be placed just above the Fabry-Perot to restrict the field 
of view from 5.7 arc minutes down to less than one arc minute as desired. The 
light passed by the 150 mm Fabry-Perot was coupled by a 3:1 ratio afocal lens 
system to a lower resolution (6A==1.5A) 50 mm Fabry-Perot which was placed in 
series with any one, of several 50 mm aperture interference filters with band- 

o 

passes typically 15-20A. The low resolution Fabry-Perot was used to suppress 
all but one of the narrow {6X=0.17A) transmission -peaks of the large Fabry- 
Perot which fell within the passband of the interference filter. The Fabry- 
Perot etalons were housed in separate gas-tight chambers and pressure-scanned 
over spectral intervals up to 6A using SF^ as the scanning gas. An automatic 
pressure difference control system maintained the tune of the two Fabry-Perot 
etalons during scans. In order to monitor changes in sky brightness and 
atmospheric transmittance, about 4% of the light incident on the spectrometer 
was directed to a reference system containing a lOOA wide filter centered near 
the wavelength being scanned. The number of photons counted by the spectro- 
meter and reference system during small equal wavelength scan intervals were 
punched on paper tape on command from an interferometric' refractometer scanned 
In unison with the Fabry-Perot etalons. This method permitted the direct 
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coniparison and addition of scans for enhancing the signal-to-nolse ratio. 

■ -in, OBSERVATIONS AND ANALYSIS 
a) Cometary Ha Spectra 

Cometary Ha 1s produced by fluorescence following solar Lg excitation 
of cometary hydrogen. Examples of scans taken for the purpose of measuring 
cometary Ha emission are shoivn in figures 1 and 2. As these scans show, 
geocoronal Ha emission (labeled component 1) was always present and in some 
cases blended with the cometary emission (labeled component 2). In order to 
isolate the purely cometary emission features, scans of geocoronal Ha emission 
were taken in directions away from the comet and away from regions of expected 
galactic emission, thus determining the shape and position of the geocoronal 
components in the comet scans. The absolute wavelength calibration was 
obtained from the position of the geocoronal Ha line, which is fluorescently 
produced by solar L3 excitation of geocoronal hydrogen. Since only the 3p 
levels are excited in this process, the center of gravity of the Ha line will 

o 

be- at 6562. 74A (Garcia and Mack 1965). In those scans were the cometary and' 
geocoronal Ha lines were well resolved (e.g., figure la) the line identified 
as cometary Ha was always shifted in wavelength from the geocoronal line by 
the amount corresponding to the comet-earth relative velocity listed in the 
ephemeris (Clark 1973) . 

Useful Ha scans were taken out to 6 arc-minutes (250,000 km) from the 
comet head; beyond this the cometary Ha was not detectable above the noise. 
Figures lb and 2b show scans taken with the field of view centered on the head 
of the comet. In these scans the cometary emission feature is stronger than 
in the sunward direction scans (figs, la and-2a) and also appears to have aii 
assymatric line profile suggesting the presence of a third emission component 
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(labeled coinponent 3) on tho red wing of the coinetary Hct line. Scans taken 
in the tailward direction out to 200,000 km (5.5 arc-minutes) from the head 
(fig. 2c) in late January show this additional component more clearly, t^e 
have tentatively identified the additional line as an H 2 O'*' emission line ' 
(Wehinger et^ al . 1974). The H 2 O emission was always red shifted v;ith respect 
to its comet rest wavelength and was only present in the scans centered on the 
head and in the tailward direction. This indicates that the was not 
expanding radially but rather was moving out in the tailward direction as would 
be expected for an ion. The measured tailward velocities of the H20'^ ions were 
in the range of 20 to 40 ktn/sec. 

In the analysis of these scans, the spectra were computer fitted using 
three emission components convolved with the instrumental profile: 1) a 

geocoronal Ha line, 2) a cometary Ha line having a width determined from the 
unblended sunward (hence presumably free) scan of 1973 December 6 and 

f 

having a position relative to geocoronal Ha determined by the comet-earth 
relative velocity, and 3) a cometary H20'^ line, red shifted from cometary Ha. 
Except for twilight scans, the effect of the solar Fraunhofer Ha absorption line 
could be neglected. _ Reasonable variations of the parameters were allowed in the 
fitting routines. In the figures, the fits obtained in this manner are shown 
superposed on the data. 

b) Absolute Ha Intensity Determination 

The determination of the comet Ha intensities (photons cm~^sec~^ ster""') 
was accomplished by two different methods. In the first method the sensitivity 
Of the spectrometer at Ha was determined by scanning a galactic emission 
nebula, NGC 1499, when it was near the zenith. The Ha surface brightness of 
a selected region of NGC 1499 was taken as 580 ± 120 Rayleighs from a previous 
comparison by Reynolds (Reynolds ^ al_^ 1973) of NGC 1499 v^ith the North 
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American nebula (NGC 7000) whose absolute Ha intensity was measured by Ishida 
and Kawajiri (1968). From several observations of the comet in a given night, 
the extinction for the large zenith angles (or large number of air masses) of 
the cometary observations was estimated, and thus comet intensities were 
obtained. This method is not very accurate due to the large and probably 
uaneliable corrections for atmospheric extinction in many of the comet scans 
taken at large zenith angles. 

The second method of obtaining absolute comet emission line intensities 
took advantage of the fact that geocoronal Ha was always present in the comet 
Ha scans. Using values of the geocoronal Ha intensities calculated from a 
theoretical model by R. R. Meier (1969, 1974), the comet Ha intensities were 
directly determined. The advantage of this method is that in each scan both 
comet Ha and geocoronal Ha undergo the same, extinction, and thus it is not 
necessary to correct explicitly for extinction. Since the success of this 
method depends on the accuracy of the geocoronal model, the model predictions 
were compared with our geocoronal intensity measurements based on the HGC' 1499 
calibration for three scans at zenith angles less than 50° where corrections 
for atmospheric extinction are reliable. On two of these scans the model 
predictions were 14% higher than the measured intensities and on the third scan 
the model prediction was 15% lower than the measured intensity. The absolute 
Ha intensities reported here were those obtained by this geocoronal Ha comparison 
method. 

Figure 3 shows the observed Ha surface brightnesses plotted as a function 
of the sun-comet distance. Since the scans were variously taken on and off the 
head and with different fields of view, the data were all normalized to yield 
the brightness averaged over a 2 arc minute field of view centered on the head 
at a distance of 1 AU from the earth. In doing the normalization, scans taken 
on the same night with different fields of view and in different positions with 
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respect to the comet head were compared for severaVseparate nights to 
estimate the Hci brightness distribution around the comet head. 


c) Hydrogen Production Rates ' - 

The observed cometary Ha brightness distribution is assumed interpretable 
in terms of a simple radial outflow model. The method of analysis used here 
to determine hydrogen production rates is similar to that in Carruthers et al . 
(1974). The hydrogen density distribution n^(r') is given by 


n^(r') = Qn/4Trr'^^ (1) 

where is the hydrogen production rate in atoms s"\ r' is the distance of 
the hydrogen from the comet nucleus, and 7^ is the average hydrogen outflow 
velocity. Then the column density H^(r) of hydrogen along a line of sight 
which has a distance r of closest approach to the nucleus, is 


Nn(r) = QH/4rv^^. (2) 


The Ha intensity along that same line of sight is given by ■ , 

I„ = gNH/4iT (3) 

where g is the number of Ha photons per second fluorescently produced by' 


each hydrogen atom in the line of sight due to absorption of solar Lyman e 
photons. This approach assumes that the hydrogen cloud is optically thin 
to Lyman g. This assumption is justified below. 

The Ha emission rate factor g is given by- 


g> 


(4) 


where g^ is the Ha emission rate factor at 1 AU and R is the sun-comet 
distance in AU. The factor is given by 

2 

where the factor 0.12 is the Ha branching ratio from the 3p state — 


31 


( 5 ) 
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is the scattering coefficient for Lyman 3, and is the solar Lyman 3 photon 
flux per unit frequency interval at 1 AU in the neighborhood of the cometary ■ 
Lyman 3 line. 

The profile of the solar L£ emission line integrated over the solar disk 
has not been measured. The L3 profile measured at one point on the sun by 
Tousey (1962) shows a self-reversed line with a peak to valley ratio of 3 to 1 . 
Geocoronal hydrogen would be excited by solar L& in the valley of this profile, 
while the comet-sun velocity during our observations Doppler-shifted the 
cometary hydrogen L0 line to a peak of the solar L3 profile. The factor 
which accounts for the geocoronal emission is 1.5 x 10"^ s“^ (Meier 1974); 
thus the factor for cometary hydrogen is taken to be 4.5 x 10“^ under 
the assumption that the Tousey profile is a reasonable approximation to the 
integrated disc profile. 

The hydrogen atom mean outflow velocity 7^ was determined from the 
observed line profiles, which had a width about two times instrumental width. ' 
The profiles cannot be fitted by assuming that the comet hydrogen atoms were 
all moving with the same speed. However, reasonable fits are obtained 
assuming that all the H atoms were moving radially outward with a speed 
distribution that was Maxwellian (Bertaux ^iL_ 1973). With this model we 
obtain a mean outflow velocity v^ of 7.8 ± 0.2 km s"'' (which corresponds to 
a temperature of 2900 ± 150°K). 

In a radial outflow model for the comet, the average brightness over 
the field of view is approximately the brightness at the center of the field 
of view If the comet head is not in the field. Analysis of the sunward scan 
in figure la (used as an example throughout this. section) gives a maximum 
hydrogen density of 1.0 ± 0.2 x 10^ at the center of the field of view 
(84,000 kin from the comet center). This corresponds to a production rate 
= 3.5 + 0.6 X 10^^ hydrogen atoms s for a comet-sun distance of 0.75 AU, 
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The effect of multiple scattering of LSwas estimated by considering a 
thermalized 2900°K hydrogen cloud with a density falling as 1/r'^ from the 
comet head and having a density of at r' = 8.4 x 10^ km.' Under 

these conditions it is readily shown that the solar L$ flux would be reduced 
by 1/e at a distance of 5.4 x 10 km from the comet head. Since the region 
within this radius would contain less than 10% of the total number of hydrogen 
atoms in the 2 arc minute field of view, we neglected multiple scattering of 
solar L3. 

The accuracy of the simple radial outflow model used here is limited by 
at least two effects (Carruthers 1974). First, solar La radiation 

pressure would tend to push the hydrogen cloud away from the sun. For the 
outflow velocity of 7.8 km s ^ found from the scan shown in figure la, one 
finds that the hydrogen atoms took about 10^ seconds to move from the comet 
head to distances corresponding to the center of the field of view. Radiation 
pressure would have slowed these atoms by only 0.07 km sec"'' , and thus would 
have had little effect on results deduced from this scan. The error in 
neglecting this effect is less than 5 % of the Ha surface brightnesses deter- 
mined from any of our scans. A second effect is the finite lifetime of the 
hydrogen cloud due to ionization by the solar wind. Carruthers et al. (1974) 
give a lifetime of about 2 x 10^ sec for the comet-sun distance for figure la, 
which means that the hydrogen atoms on the average passed well beyond the 
region explored in this study before the effects of ionization v/ere significant 
The error in neglecting this effect on the surface brightness is less than 5% 

■ for any of our scans. 

Figure 4 shows the hydrogen atom production rates plotted versus the 
comet-sun distance. Scans centered off the comet head were used to determine 
the absolute production rate in the manner indicated above. A direct deter- 
mination of the production rate from fiead-centered scans would require a model 
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for hydrogen densities near the nucleus. To avoid this problem, only 
relative production .rates were determined from' the ' head-centered scans/and 
these in turn were normalized to the production rates determined from scans 
off the head on those days in which scans both on and off the head were made.' 

d) Deuteri urn /Hydrogen Ratio 

In an attempt to measure the D/H ratio, several scans were extended to 

o 

include the region 1.8 A to the blue of Ha which is the expected position of 
the Da emission line. Figure 5 shows preperihelion and postperihelion scans. 
Before perihelion the comet-sun velocity Doppler shifted the cometary hydrogen 
L3 line to the red peak (assuming the Tousey profile) of the solar' L6 line and 
the cometary deuterium L3 line to the blue peak of the solar L8 line. After 
perihelion the solar spectrum was red shifted relative to the comet, placing 
the cometary hydrogen LB line on the blue peak of the solar LB line and the 
cometary deuterium LB line completely off the solar LB Tine. Thus postperihelion 
observations should show no Da feature. In preperihelion observations cometary 
deuterium and hydrogen would have approximately the same g factor and there- 
fore a direct D/H ratio could be measured if the Da line were detected. 

Scans for Da were taken to within 10 days of perihelion. The spectrum is 
complicated by two terrestrial H 2 O absorption features, a weak cometary H 20 '^ 
line and possibly other weak airglow or comatary features. In the near- 
perihelion period the comet was visible only in the morning twilight sky where 
terrestrial features were enhanced. Ten days before perihelion the morning 
twilight dominated even the cometary Ha emission and measurements were halted. 
After perihelion, observations in the Da region of the comet spectrum were made 
for comparison. Observations of the morning twilight were also made to remove 
terrestrial features. An analysis of all Da scans did not show any significant 
additional feature on the preperihelion scans that could be attributed to 
cometary Da emission. Based on the noise in the Da region we are able to put 
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an upper limit of 1 part in 100 for the D/H ratio in the .cometary atomic 
hydrogen cloud. 

e) [01] X630Q A Emission 

By changing the interference filter and retuning the etalons, observations 

' o ' ■ . ■ 

could also be made at 6300 A in order to observe [01] emission from the comet. 

By comparing filter efficiencies, etalon transmission curves and photomultiplier 
quantum efficiencies, the instrument efficiency at 6300 A relative to the 
efficiency at Ha could be determined, permitting absolute intensity measurements 
Atmospheric extinction estimates were made by comparing several scans made on 
the same day, as in the first method described above for Ha. 

o 

All 6300 A scans, examples of which are shown in figure 6, contain the 
[01] airglow feature at 6300.23 A (Chamberlain 1961). The measured airglow 
surface brightness indicates that the calibration was reasonable; however it is 
not sufficiently well known to use for calibration. Scans on and near the comet 
head show an additional line shifted from the airglow line by the Doppler 
shift corresponding to the comet-earth relative velocity, and is identified as 
cometary [01] X6300 emission. This cometary [01] emission presumably is 
produced by molecular dissociations which leave the oxygen atoms in the 
state. Collisional de-excitation occurs only in a very small region near the 
center of the comet, and both fluorescence and collisional excitation must 
be small (Arpigny 1965, Blermann and Trefftz 1964). 

_ The cometary line was consistently slightly wider than the airglow line, 
but the uncertainty in the derived width is large because the width was 
approximately one half the instrumental width, which was chosen to be optimum 
for the hydrogen measurements. The measured full width at half maximum of the 
cometary [01] profile was 2.7 (+1.0. -1.3) km (corresponding to an effective 
temperature of 2400 [+2200, -1300]°K). These velocities are insufficient to 
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carry the excited atoms outside the field of view before radiating. There- 
fore, the absolute intensity measurements of [01] X6300 give directly the 
production rate of oxygen atoms in the state. The surface brightness of 

200 ± 100 Rayleighs measured on 1973 December 8 gives a oxygen production 

28 “11 

rate of 2 ± 1 x 10 s . The D oxygen production rates have been plotted 
as- a function of comet-sun distance in figure 7. 

IV. DISCUSSION 

The values v/e have obtained for the hydrogen production rate at various 
heliocentric distances, shown in figure 4, are generally in good agreement 
with the results obtained by Carruthers et al_^ (1974) using Lyman alpha 
imagery of Comet Kohoutek. , At present, the major source- of. uncertainty in 
our results for hydrogen is probably the lack of reliable solar Lyman beta 
profiles,, -integrated over the solar disk, for different levels of solar 
activity. Presumably, such observations will eventually be carried out. 

Additional analysis of the red shifted K20'^ emission line detected in 
our Ha scans is being carried out in order to obtain an estimate of the H 2 O'*' 
production rate in Comet Kohoutek. It should then also be possible to 
obtain an estimate of the H 2 O production rate, using the relevant cross- 
sections for photodissociation and photoionization of H 2 O. 

Since a number of different molecules, including H 2 O, could be sources 
of the hydrogen and oxygen atoms which produced the emissions we observed, 
we cannot at present use our observations to determine the production rate 
of oxygen or any of the possible parent molecules. However, the b oxygen 
production rates at least set a lower limit to the total oxygen production 
rate, and thus our results are consistent with the observations of Feldman et 
al- (1974) who found a production rate of 1.4 x 10^^ s“^ for oxygen. These 
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rates indicate that about 14% of the oxygen atoms were produced directly into 
singlet states. 

In conclusion, we feel that ground-based Fabry-Perot observations of 
comets can provide a powerful tool for the study of comets, especially when 
carried out in close coordination with other types of observing techniques. 

We would like to thank F. Barmore and B. Donn for their assistance and 
advice. We are especially grateful to R. R. Meier for calculating the geo- 
coronal Ha intensities used in our cometary Ha intensity determinations. 
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FIGURE CAPTIONS 

Fig. 1 Spectral scans near Ha taken on 1973 December 6 with a circular 2 
arc-minute field of view. Scan (a) was taken at 1250 UT with the 
center of the field 1.5 arc-minutes (84,000 km) sunward of the comet 
head. The geocoronal Ha line (component 1) is at 6562.74 A, and the 
'' cometary Ha line (component 2) is blue shifted 0.63 A corresponding 
to the comet-earth relative velocity. The continuous line through 
the data is generated by a continuum and a two-Gaussian fit, the 
components of which are shown below the data. Scan (b) was taken at 
1235 UT centered on the head. A third emission feature is evident, 
and the fit shown includes a third Gaussian component (component 3) 

o 

0.3 A redward of cometary Ha. Scan (c) was taken at 1300 UT with the 
field of view centered 1.5 arc-minutes tailward of the comet head. 

The' fit again includes a third component red-shifted by 0.3 A from 
cometary Ha, and identified as HgO'*’.. 

Fig. 2 Spectral scans near Ha taken 1974 January 26 with a circular 6 arc- 
minute field of view. Scan (a) is the sum of scans taken at 0254 UT 
and 0428 UT and centered 5.5 arc-minutes (2.1 x 10^ km) sunward of 
the comet head. The cometary Ha line (component 1) is red-shifted 
from the geocoronal Ha line (component 2) at 6562.74 A due to the 
comet-earth relative velocity. Scan (b), taken at 0350 UT centered 
on the comet head shows the feature (component 3). Scan (c) is 
the sum of scans taken at 0307 UT and 0416 UT with the field centered. 
5.5 arc minutes tailward of the comet head. The third component is 
again evident. 
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Fig. 3 Brightness of cometary Ha emission. The ordinate gives the brightness 
averaged over a 2 arc-minute head-centered field at a distance of 1 AU 
from the earth as a function of the comet-sun distance. Preperihelion 
data are shown as open circles and postperihel ion data are shown as 
closed circles. 

Fig. 4 Hydrogen atom production rate as a function of the comet-sun distance. 

Fig. 5 Comparison of preperihelion and postperihelion scans near Ha, including 
the region where Da emission would occur. Positions marked are (from 
left to right): unidentified terrestrial airglow emission, cometary 
H 2 O'*' emission, cometary Da emission, terrestrial H 2 O absorption, 
cometary Hci emission, and geocoronal Ha emission. 

O 

Fig. 6 Spectral scans taken near [01] X6300 A on 1974 January 24 with a 

circular 40 arc-second field of view. Scan (a) was taken at 0230 UT 
centered on the comet head, and shows the cometary [01] X6300 A 
emission red-shifted due to the comet-earth Doppler' shift by 0.45 A 
from the [01] airglow feature at 6300.23 A. Scan (b) shows the same 
spectral region scanned at 0240 UT with the field of view centered 40 
arc-seconds (24,000 km) off the comet head. Scan (c) is the sum of 
two scans taken at 0205 UT and 0250 UT with the field centered 80 arc- 
seconds off the comet head. 

Fig. 7 Production rate of D oxygen atoms as a function of the comet-sun 
distance. 
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